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Abstract. Broad-band (gamma to radio) variations of the flux 
density were observed in the first half of 1992 in the luminous 
high redshift (z = 2.172) quasar S5 0836+710. VLBI moni- 
toring observations during 1993 - 1996 performed at 86 GHz, 
22 GHz, 15 GHz, and 8 GHz show the ejection of a new jet 
component, which most probably is directly related to a quasi 
simultaneous gamma-, X-ray, optical flaring activity which was 
observed in February 1992. During the period 1992 - 1993 the 
flaring propagated through the radio spectrum. From several 
quasi-simultaneous radio spectra taken during this phase of ac- 
tivity, we determine the time evolution of the spectral turnover 
of the radio spectrum in the S m -v m diagram. The data indi- 
cate a correlation of the jet activity with the variability of the 
broad-band electromagnetic spectrum of the source. The obser- 
vational findings are discussed in the framework of relativistic 
shock models. 
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1. Introduction 

In early 1992, the ultraluminous S5 quasar 0836+710(4C71.07, 
z — 2.172) underwent a prominent optical outburst (von Linde 
et al. 1993), which lat er led to enhanced variabili ty in the mm- 
and cm- radio bands ( [Marscher and Bloom 1994 ). At the time 
of the optical flaring, the source also was in a bright state in 
the gamma-regime (EGRET: Fichtel et al. 1994) and at soft X- 
rays (ROSAT: [Brunner et al. 1994 ). Motivated by the increased 
activity of the source, we started a high frequency VLBI mon- 
itoring program in 1993 in order to investigate the possible 
relation between broad-band changes in the flux density and 
structural variations in the jet on sub-milliarcsecond scales. 
The ejection of new VLBI jet components after major out- 
bursts in the flux density has been observed also in a few other 
sources (e. g. component ejection after optical flaring: 3C 345 



Babadzhanyants and Belokon 1986, 3C 273 Krichbaum et al. 1990a; 



ejection after radio flares: BL Lac Mutel et al. 1990, PKS 0420-014, 
Wagner et al. 1995a) and seems to be a quite common phe- 
nomenon in compact flat spectrum radio sources (blazars). 

With the new data obtained for 0836+710, it appears very 
likely that the relation between outbursts in the flux density ob- 
served at high frequencies and jet component ejection is not 
limited to the synchrotron regime (e.g. optical-to-radio), but 
most probably covers a much wider spectral range (gamma-to- 
radio), possibly the full electromagnetic spectrum. Similar cor- 
relations between gamma-activity and jet activity were recently 
found also in a number of other sources, e.g. in 3C 273 (Krich- 
baum et al. 1996), 3C 279 (Wehrle et al. 1994), 0528+134 
and 3C454.3 (Krichbaum et al. 1995), where the ejection of 
jet components could be related to gamma-flares, or enhanced 
levels of gamma-activity. However, in most of the published 
cases the frequency coverage is not as broad as presented here 
for 0836+710. 

The quasar 0836+710 is a member of a complete flux den- 
sity limited sample of 13 flat-spectrum radio sources (Eckart 
et al. 1986 & 1987) compiled from the 6 cm S5 survey (Kiihr 
et al. 1981). The 13 sources of the sample are located in the 
northern polar cap (S > 70°) and are regularly studied with 
VLBI at frequencies ranging from 0.3 to 22 GHz (cf. Witzel et 
al. 1988 and references therein). The VLBI data for 0836+710 
now span a time range of more than 15 years (1983-1996). 

The high dynamic range VLBI images of 0836+710 reveal 
a complex and wiggled one-sided core-jet structure (Eckart 
et al. 1986, [Eckart et al. 1987] ). World-array VLBI images at 
90 cm and 18cm showthejet extending overmorethan 150milliarc- 
seconds (mas) in direction to the outer arcsecond lobe seen with 
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MERLIN and the VLA ( |Hummel et al. 1992[ ). Mainly from the 
VLBI monitoring at 6 cm (Krichbaum et al. 1990b hereafter 
K90b, Otterbein 1996 hereafter 096), a very complex motion 
pattern is revealed, with jet components, moving at velocities 
ranging from subluminal (/3 ap p = 0.5; /3 = v/c) to super- 
luminal (/3 a pp = 10) (assuming Ho = lOOkm/s/Mpc and 
qo = 0.5, throughout this paper). No systematic correlation be- 
tween the speed of the jet and the separation of the components 
from the core seems to be present. The jet shows a noticeably 
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Table 1. Summary of experiments 



date 



freq. 
[GHz] 



observing mode, 
bandwidth [MHz] 



array 



1993, Apr. 7 
1993, Sep. 16 
1995, Jan. 20 
1995, Aug. 24 
1995, Aug. 24 



86.2 MKIII(A), 112 

22.2 MKIII(B),56 

15.3 VLBA 64-8-1, 32 
22.2 VLBA 128-8-1,64 

8.4 VLBA 128-8-1,64 



global (4) 1 
global (11) 
VLBA (10) 2 
VLBA(IO) 3 
VLBA (10) 3 



Notes: 

1 : stations: Effelsberg, Pico Veleta, Onsala, Haystack 
2 : snapshot-type observations 

3 : 8/22 GHz simultaneous dual-frequency observations 



structure, with lateral displacements of its ridge line (kinks) and 
oscillations of its transverse width (K90b). 

At cm-wavelengths, 0836+710 is strongly polarized (p ~ 
9 % at 5 GHz). Polarization VLBI images show a highly po- 
larized jet with the magnetic fiel d roughly following the bent 
ridge line (Cawthorne et al. 1993). This suggests that the jet of 
0836+710 is also highly magnetized. 

In this paper, we focus on the structural changes observed 
in the sub-milliarcsecond regions of the inner jet of 0836+710 
during and after the optical outburst of 1992. A more detailed 
discussion of the overall morphology and the spectral and kine- 
matic properties of the milliarcsecond-jet will be given else- 
where. 
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2. Observations and data reduction 



Fig. 1. Overall spectrum (jaF„) of 0836+710 as compiled from 
presented and published data. 



86 GHz we included corrections for atmospheric absorption (opac- 
ity) in the amplitude calibration. 

The imaging was performed using the CalTech VLBI pack- 
age and the difference mapping software DIFMAP (Pearson 
1991, Shepherd et al. 1994). In parallel to the mapping, we 
made Gaussian model fits to the visibility amplitudes and phases, 
in order to parameterize the multi-component structures seen 
in the maps. The measurement errors of the parameters of the 
Gaussian components were determined formally from least square 
fits and, in addition, from their variations on fits to slightly dif- 
ferent calibrated and edited data sets. 



0836+710 was observed with VLBI between 1993 and 1995 
at frequencies ranging from 8 GHz to 86 GHz. The observa- 
tions in August 1995 were made in dual-frequency mode (22 
and 8 GHz). At 86 GHz and 15GHz[]the data for 0836+710 
were obtained during VLBI campaigns which also included 
other sources (86 GHz: Schalinski et al. 1994; 15 GHz: Pat- 
naik et al. 1996). Owing to the limited uw-coverage in the 86 
and 15 GHz data sets, the dynamic range of the corresponding 
maps is lower than that of the high dynamic range images at 22 
and 8 GHz, with their optimal uv-coverage. The details of the 
observations are summarized in Table [l]. Column 1 gives the 
observing date, column 2 the frequency, column 3 the observ- 
ing mode and bandwidth, and column 4 the VLBI array with 
the number of participating stations in brackets. 

After the correlation (VLBA data in Socorro, MK III data 
in Bonn), the data were fringe-fitted and calibrated in the stan- 
dard manner. AIPS was used for reducing the VLBA data. The 
fringe fitting of the MK III data was done using the Haystack 
software (Rogers et al. 1983), and the baseline oriented global 
fringe fitting routines available at the MPIfR (Alef & Pore as 
1986, Alef 1989). The amplitude calibration was carried out 
in the standard way, using system temperature measurements 
and elevation-dependent antenna gain curves. At 22 GHz and 

1 The 15 GHz data were kindly provided by A. Patnaik 



3. Results 

3.1. Broad-band Variability 

0836+710, which in the past showed only mild variations of its 
radio flux density, recently changed into a more active phase 
of variability. Between 1990-1993 the source exhibited strong 
variability, which seems to be correlated over a wide range of 
the electromagnetic spectrum. In Figure [T], we show the total 
luminosity spectrum; and in Figure g, we show the flux den- 
sity measurements for the interval 1990.5-1993.2. The data 
are taken from the literature and include, from the top to the 
bottom, EGRET observations in the gamma-band at 100 MeV 
(Mukherjee et al. 1997), ROSAT observations in the X-ray band 
at IkeV (Brunner et al. 1994), optical R-band observations 
(von Linde et al. 1993, Schramm et al. 1994 and Heidelberg 
telescope), 230 GHz and 90 GHz flux density measurements 
from the IRAM 30 m millimeter telescope at Pico Veleta and 
from the JCMT at Mauna Kea (Steppe et al. 1993, Bloom et 
al. 1994), and measurements at 43, 37 and 22 GHz from the 
Metsahovi flux density monitoring program (E. Valtaoja, priv. 
comm.) and from the 100 m radio telescope at Effelsberg. 

From Figure |[ it is clear that the sparse and irregular time 
sampling, particularly in the X-ray and gamma-bands, requires 
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Fig. 2. Broad-band flux density variability in 0836+7 10 for the 
period 1990.5 - 1993.2. From top to bottom gamma-ray, X-ray, 
optical, millimeter- and cm-radio data are shown (for origin of 
data see text). The lines connect adjacent data points, arrows in 
the top panel indicate upper limits. The data suggest a corre- 
lated gamma- to optical outburst first peaking in 1992.13, and 
propagating later towards longer (mm-/cm-) wavelengths. 



a careful interpretation. Obviously, a cross correlation analy- 
sis, or an accurate determination of the time lags between the 
different bands is not possible. The data nevertheless indicate 
a tight correlation between the gamma-/X-ray and optical flux 
density levels. Close to the optical flaring activity observed in 
1992.1-1992.2, the gamma- and X-ray flux density was also 
high. The EGRET measurements show variability by a fac- 
tor of 2-3 on a time scale of < 0.5 yrs, and low flux den- 
sities shortly before and after the time of the optical flaring. 
The X-ray flux is also variable. About 0.5 years after the op- 
tical outburst, it dropped by a factor of 2. Given the sparse 
sampling, the data are consistent with a (at least nearly) si- 
multaneous (At < 0.1 yrs) gamma-ray/optical flaring. This is 



consistent with nearly simultaneous gamma-/optical variations 
observed in a few other objects, eg. in 1406—076 (Wagner et 
al. 1995b). From the measurements in the mm- and cm-bands 
it appears that the optical outburst propagates towards longer 
wavelengths, with the maximum in the radio appearing after 
0.1 — 0.3 yrs first at shorter (mm-), and after 0.3 — 0.5 yrs later 
at longer (cm-) wavelengths. 



3.2. Results from the VLBI monitoring 

3.2.1. Overall morphology and component identification 

In Figure [|, we show two 22 GHz images of 0836+710 ob- 
tained in 1993.71 and 1995.65. As in the maps obtained earlier 
at longer wavelengths, the source shows a one sided core/jet 
structure with several compact components. From the maps and 
the model fits, we identified several distinct jet components B3, 
B2, B, and CI (B and CI are blended) located at separations 
from the core of up to r ~ 4 mas. The labels in Figure || reflect 
this identification scheme, in which we implicitly assumed that 
the brightest component located at the map center is the sta- 
tionary core of the jet. Beyond a separation of 4 mas from the 
core, the jet emission becomes weaker and more extended, so 
that it is resolved by the interferometer beam. Images obtained 
by applying a taper to the 22 GHz uv— data, however, show 
that the emission from the jet extends to 15 mas from the core, 
in good agreement with the outer structure seen previously at 
longer wavelengths (K90b, Hummel et al. 1992). 

In the following, we use a component annotation similar 
to the one introduced by K90b for the 5 GHz data. The fea- 
ture seen at a separation of r ~ 3 mas from the core can be 
identified as a blend between a fast moving older 5 GHz com- 
ponent B (/9 a pp = 9.6 ± 0.2; /3 = v/c) and a quasi-stationary 
component CI (/3 app < 2). Near component CI, the jet di- 
rection changes. This displacement is seen more clearly in the 
'tapered' 22 GHz images, and also in the maps at 15 and 8 
GHz. We note that similar displacements also appear at larger 
separations from the core, in particular near component D at 
r ~ 10 mas in the 5 GHz maps (K90b), and near component 
G2 at r ~ 30 mas (Hummel et al. 1992) at 1.6 GHz. With the 
new 22 GHz maps, it appears very likely that the lateral dis- 
placements (kinks) of the jet axis are systematic and reflect a 
general jet property in 0836+710. 

At smaller separations from the core, two new components 
B2 & B3 are seen (Figure |]). Linear back-extrapolation of the 
motion of B2 allows its identification with a component pre- 
viously visible in a 5 GHz map of 1990 (096). Adopting its 
constant apparent speed of /3 app — 8.3 ± 2, the component B2 
was ejected before 1986.8 and therefore cannot be related to 
the flare observed in 1992. 

In the 22 GHz map of 1993.71, the core appears elongated, 
indicating the emergence of a new component (see Fig. || left). 
From Gaussian model fits we determined its relative separation 
from the core of r ~ 0.3 mas. We label this component B3. 
The image of 1995 (Fig. |, right) shows B3 at r ~ 0.8 mas, 
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Fig.3. 22 GHz maps of the inner jet of 0836+710 from 1993.71 (left) and 1995.65 (right). Contour levels are -0.25, 0.25, 0.5, 
1, 2, 5, 10, 20, 30, 50, 70, 90 % of the peak flux density of 0.35 Jy/beam (left), and 0.9 Jy/beam (right). Both maps are restored 
with a circular beam of 0.45 mas size. Labels denote the component identification. 
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Table 2. Relative separation from the core of the new compo- 
nent B3 



epoch 


freq. 


R 


PA 




[GHz] 


[mas] 


n 


1993.27 


86 


0.15 ±0.05 


220 ± 20 


1993.71 


22 


0.34 ± 0.04 


231 ±3 


1995.05 


15 


0.72 ±0.15 


227 ± 5 


1995.65 


8 


0.69 ±0.15 


229 ± 5 


1995.65 


22 


0.76 ± 0.06 


230 ± 3 


1996.34 1 


22 


1.05 ±0.10 





Notes: 1 : A. Marscher, priv. comm. 



Fig. 4. 86 GHz map of 0836+710 of 1993.27. Contour levels 
are -1, 1, 2, 5, 10, 20, 30, 50, 70, 90 % of the peak flux density 
of 0.41 Jy/beam. The observing beam size is 0.30 x 0.06 mas, 
P.A. = 75°. 



thus indicating superluminal motion with a speed similar to its 
predecessors B2 and B (see next paragraph). 

The ejection of a new component and its subsequent sep- 
aration from the core is further supported by mm- VLBI ob- 
servation of the source. A high angular resolution map ob- 
tained at 86 GHz in 1993.27 (Fig. ||), 5 months earlier than 
our first 22 GHz observation, showed a secondary component 
next to the unresolved VLBI core (size < 0.1 mas). Owing to 
the small number of participating VLBI stations at this high 
observing frequency and the resulting limitations of the uv- 
coverage, the 86 GHz map shows only the most prominent fea- 
tures of a probably more complex brightness distribution. The 
secondary component is located at r ~ 0.15 mas. Since its po- 



sition agrees well with the extrapolated position of B3 (from 
22 GHz data) at that time, we identify this component with B3. 



3.2.2. Motion of component B3 

In Table ||, we summarize the relative separations of B3 from 
the core for all available observing epochs. We further added a 
data point from a 22 GHz map of 1996.34, which was kindly 
provided by A. Marscher (priv. comm.). 

In synchrotron self-absorbed VLBI jets (e.g. Blandford and 
Konigl 1979), frequency dependent shifts of the relative posi- 
tions of the VLBI components, caused by optical depth effects, 
are expected. The comparison of the relative component posi- 
tions between the nearly simultaneous 8, 15 and 22 GHz im- 
ages, however, indicate that the position shifts in 0836+710 are 
relatively small (^ 0.1 mas). In the following, we neglect such 
shifts due to opacity and plot in Figure || the relative separation 
from the core of B3 versus time. The data points line up very 
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Fig. 5. Relative separation of component B3 from the core 
plotted versus time. Symbols denote different frequencies: di- 
amonds (8 GHz), triangles (15 GHz), circles (22 GHz), and 
squares (86 GHz). The solid line represents a linear fit with 
slope /it = 0.26 mas/yr, corresponding to superluminal mo- 
tion with /3 app = 10.8. Linear back-extrapolation (dashed line) 
yields component ejection near 1992.65. The arrow indicates 
the onset of the optical flare (1992.13). 



nicely, thus confirming indirectly the validity of our assumption 
of negligible opacity effects. Furthermore, we find no evidence 
for non-linear motion, e.g. acceleration or motion along curved 
paths (see also the nearly constant position angle of B3 in table 

A linear fit to the data yields an apparent angular separa- 
tion rate of /i = 0.26 ± 0.03 mas/yr corresponding to an ap- 
parent velocity of /3 a p P = 10.8 ± 1.3. This velocity is close 
to the speeds seen for the older components B & B2 (K90b, 
096), and confirms the previous finding of relatively high ve- 
locities at small separations from the core (note: components 
CI and C2 move substantially slower, cf. K90b). Linear back- 
extrapolation of the motion of B3 yields the time of zero sep- 
aration of < C j = 1992.65^Q'}y. This is close to the time of the 
strong optical flare observed in 1992. 13 (von Linde et al. 1993, 
see Figure |J). Within possible extrapolation errors and due to 
the limited and partially incomplete time coverage of the opti- 
cal flaring, it therefore appears that B3 was ejected during or 
shortly after (0 ^ At < 0.67 yrs) the optical outburst. We 
study the correlation between the optical flaring and the ejec- 
tion time of B3 in more detail in section 4.2 



3.3. Evolution of the total radio spectrum 

Figure || shows the propagation of the outburst through the 
spectrum. Such a behavior is readily described within the stan- 
dard shock-in-jet models (e.g. Marscher & Gear 1985), in which 
a relativistic shock propagates along a synchrotron self-absorbed 
jet, causing (at least qualitatively) the observed frequency de- 
pendence. 
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Fig. 6. Spectrum of 0836+710 at epoch 6 (1992.75). In the up- 
per picture the total spectrum is plotted. The dashed line shows 
the quiescent spectrum, i.e. the quiescent part determined from 
an observations prior to the optical flare in 1992. In the lower 
panel the residual variable part of the spectrum (quiescent spec- 
trum subtracted) and the fit used to determine the spectral turn 
over is shown. 



From the data shown in Figure g and additional data from 
measurements made with the MPIfR 100 m radio telescope at 
longer cm-wavelengths, we are able to follow the evolution of 
the radio spectrum and to determine the path of a local spectral 
excess, which propagated from higher to lower frequencies. 

From our radio data, we first constructed quasi-simultaneous 
radio spectra covering the 1.4 — 230 GHz range. Good spectral 
coverage is achieved in the period from 1992.05 to 1993.00. To 
overcome the lack of frequency coverage at individual epochs, 
we binned the available data into 7 intervals, each of 50 days 
in duration. Since the source did not show strong variations 
within the binning-timescale, we are confident that the averag- 
ing of the data within each bin is not causing large errors for our 
later derived turnover points. Figure ^ (top) gives an example 
of such a radio spectrum. 
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Fig. 7. Path of the spectral evolution of the turnover point in the 
S m — i^m-diagram. The labels 1-7 denote the individual spectral 
epochs (see text). The straight line shows a power law fit (S m oc 
^ 58 ) to the data. 



ation, the Compton, synchrotron, and adiabatic phase. During 
the first two stages, the turnover flux rises with decreasing fre- 
quency, then remains nearly constant. The phase of adiabatic 
losses is characterized by a decrease of the flux density with 
decreasing turnover frequency. We therefore identify the spec- 
tral evolution seen in Figure [7| with the adiabatic-loss stage. A 
power law fit (the straight line in Figure 0) to the path of the 
turnover point yields S m oc i^ 5S - 

We note that in Figure a slightly steeper slope is ob- 
tained, if the first data point (point no. 1) is not included in 
the fit. Within the afore mentioned model, it is possible that the 
spectral evolution from point no. 1 to no. 2 reflects the tran- 
sition from the synchrotron to the adiabatic cooling phase. If 
this is the case, the synchrotron loss dominated phase would 
have ceased t aync < 30 days after the onset of the optical flare. 
Obviously observations much closer to the time of the outburst 
would have been required to study the synchrotron phase. 

4. Discussion 



The long-term flux density monitoring of 0836+710 at cm- 
wavelengths shows longer periods of low or only moderately 
slow (timescales of ^years) variability. One such period had 
occurred prior to the time of the optical activity in 1992. We 
used the spectra from this period to determine the shape of the 
quiescent spectrum (non-variable part). Throughout the radio 
band considered (~ 1-100 GHz) it can be well represented by 
a simple power law {S^ oc v~ a - 7 ), which we plot as dashed line 
in Figure || (top). 

In order to determine the position {v m , S m ) of the turnover 
point of the variable part of the radio spectrum, we now sub- 
tracted the same quiescent spectrum from our binned spec- 
tra. The fit of a homogeneous synchrotron component spec- 
trum (with spectral indices athick = 2.5 and athin = —0.7) 
to each of the residual spectra was rather unsatisfactory, indi- 
cating that the variable part of the radio spectrum cannot be 
represented well enough by a single and homogeneous syn- 
chrotron component. A much better fit was achieved when we 
used cttinck = 1-0 and a t hm = —0.7. In the following we use 
the results from these fits. Figure || (bottom) shows an exam- 
ple of such a fit for one of the spectra (epoch 6). From each of 
the fits to the 7 residual spectra, we determined the position of 
the spectral turnover. In Figure [7] the turnover flux density S m 
is plotted versus turnover frequency z/ m as a function of time. 
The numbered labels at each position correspond to the follow- 
ing epochs: 1 - 1992.07 (TJD 8650), 2 - 1992.21 (TJD 8700), 
3 - 1992.35 (TJD 8750), 4 - 1992.48 (TJD 8800), 5 - 1992.62 
(TJD 8850), 6 - 1992.76 (TJD 8900), 7 - 1992.90 (TJD 8950). 

The measured path of the spectral turnover seems to follow 
a linear trend (in the log(S m ) / log(u m ) space), which can be in- 
terpreted within models invoking a shock moving along a rela- 
tivistic synchrotron self-absorbed jet. In the model of Marscher 
& Gear (1985) and Marscher (1990), the spectral evolution is 
described by three subsequent phases of energy loss and radi- 



4.1. The total spectrum and the VLBI jet 

The S m — v m dependence determined in the previous section 
can be used for studying the relation between the spectral and 
kinematic properties of the relativistic jet in 0836+710. De- 
spite some evidence for emission from an inhomogeneous syn- 
chrotron component (a t hick = 1.0, see previous section), we 
postulate for the sake of simplicity, that the observed spec- 
tral variations are produced by a single, compact relativistic 
shock dominating the source radio emission, and associated 
with the core of the jet or an emitting region moving inside the 
jet within ^ 0.1 mas distance from the core (a distance com- 
parable to the back-extrapolated separation of B3 at the last 
spectral epoch, 1992.9). We use the formalism developed by 
Marscher & Gear (1985), and account for possible variations 
of the Doppler factor in the emitting region (Marscher 1990). 
Based on the synchrotron spectral index a = —0.7 used for the 
spectral fitting, a power law electron energy distribution with 
s = 2.4 (s — 1 — 2a) is assumed. The magnetic field distri- 
bution along the jet axis r is defined by an exponent a, so that 
B oc r _a . The Doppler factor is 6 oc r b . The spectral evo- 
lution of the shock is then described by the variations of the 
turnover point: S m oc i/£ and z/ m oc r e . The exponents p and e 
depend on the dominating type of the energy losses (Marscher 
1990, Marscher et al. 1991, Lobanov & Zensus 1998). Follow- 
ing Marscher (1990), we obtain the following expressions for 
p and e during the adiabatic-loss stage: 

(19 - 4s) - 3a{2s + 3) + 36(3s + 7) 
P ~ 2(2.s + l) + 3(a-&)(s + 2) ' () 

2(2s+l) + 3( a -6)(s + 2) 

e = WTl) ■ (2) 

We now consider two cases: a) a jet with transversal (a = 1) 
magnetic field, and b) a jet with longitudinal (a = 2) magnetic 
field. With the measured p = 0.58, changes of the Doppler 
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factor and turnover frequency along the jet can be determined 
from tfty and (g). 

For a = 1, we obtain b = 0.0 and e = —1.3. The dis- 
tance travelled by the emitting region between 1992.07 and 
1992.90 (the 1st and the 7th spectral epochs) should increase 
by a factor of r 7 /ri w 2.1. To assess the je t kinematics, we 
use /3 app = 10.8 measured in Section 3.2.2 for the jet com- 
ponent B3 at small separations from the core. For now, we 
assume that the jet bulk Lorentz factor is 7j = 10.9, with 
the resulting Doppler factor <5j = 11.4 and jet viewing angle 
8j = 5°. Using the time separation between the two afore- 
mentioned spectral epochs, At = 0.83 yrs, one can then esti- 
mate the distance between the emitting region and the jet apex, 
n = (3j8At/(l + z) = 9.5 ± 0.8pc. 

The case with a = 2 results in b = 0.4, implying that the 
emitting region was accelerating. The corresponding turnover 
frequency evolution is described by e = —1.7. The resulting 
Doppler factor must increase between t\ and t-j by 87/81 = 
(f-m 7 /Vm i) b ^ e ~ 1-3, and the corresponding distance must be- 
come larger by a factor of « 1.8. For a jet with variable Doppler 
factor, the distance travelled by the shock between two epochs 
t\ and <2 is (Lobanov & Zensus 1998): 
rt * (3(t)dt 



An,2 = (l + z) 



1 



(3) 



f3(t)cos0{t) ' 

For the Doppler factor 8 oc r b , equation gives for the dis 
tance to the jet apex at the epoch t\\ 



l + z 
SicAt 



1 



dr 



1/(6-1) 



(4) 



where c is the light speed, A< = £2 — ii,andr n = (f m 2/i / mi) 1 ' e 
The exact form of j(r) is essentially unknown. In many cases, 
it can be assumed to be constant (so that the Doppler factor 
variations are entirely due to a curved path followed by the 
emitting region). For a straight jet, the j(r) can be derived 
from the Doppler factor variations. For 0836+710, both ap- 
proaches yield similar result: r\ ~ 30 pc (here we postulate 
87 = 8j). One can argue that the location r-y (the beginning of 
the adiabatic-loss stage) should correspond to the VLBI core, 
since the core emission is likely to be dominated by Comp- 
ton and synchrotron energy losses (Unwin et al. 1994, 1997), 
whereas the moving jet components are successfully modelled 
as adiabatically expanding spherical plasmons embedded in the 
jet (Zensus et al. 1995). In view of this argument, the case of a 
non-accelerating jet with r\ w 10 pc appears to be more likely 
to explain the observed spectral behavior in 0836+710. 



4.2. Travel times and radio-to-gamma correlations 

If the jet core is located at 7*1, one can calculate the travel time 
between the jet nozzle, r , and the core. For b — (8 =const), 
the calculation is trivial, and it gives (in the observer's frame) 
Ai ,i = (1 + z)ri7" 1 (5- 1 /3- 1 = 0.74 ± 0.07yrs. For an ac- 
celerating jet, the travel time can be calculated from (^J), and 
it depends on the (unknown) Doppler factor 8q at the jet noz- 
zle. If the nozzle is formed by a pressure drop in a relativistic 



outflow, the resulting Lorentz factor is 70 ~ 1.2 (Marscher 
1980). The corresponding Doppler factor is then 8q « 1.8, for 
9q = 9j = 5°. The resulting travel time is Ato,i = 1.6yrs. For 
larger 70, the travel time increases and may become as large as 
2.8 yrs, if 70 ~ 71. 

On the basis of the arguments presented in the previous 
section, we take A£q,i = 0.74 yrs as a better estimate of the 



travel time between the nozzle and the core. In section 3.2.2, 
we have found t c} — 1992.65 for the ejection time of the super- 
luminal feature B3. Together with the derived nozzle-to-core 
travel time, this implies that the plasma condensation responsi- 
ble for the emission of B3 had travelled through the nozzle at 
to = 1991.91 ± 0.21. The latter epoch is very close to the time 
when the optical flare was observed in the source. Although the 
coincidence may be only fortuitous, it is rather striking, and we 
would like to emphasize that such a situation is indeed possible. 
For instance, if a flare occurs when a dense plasma conden- 
sation travels through the jet sonic point, the optical depth in 
the radio bands can exceed unity, and only a high-energy flare 
(optical, X-rays, or gamma-rays) is observed. An increase of 
the radio emission and the associated ejection of a new VLBI 
feature are observed later, after the condensation has travelled 
outside the r > 1 region. 

The jet geometry and travel times above are determined for 
7j w 7 m ; n . If we suppose that the gamma-to-optical flaring ac- 
tivity and the ejection of VLBI component B3 in 1992 had a 
common origin, the jet geometry can be uniquely determined. 
We assume that the maximum of the flare occurs at ifl aro = 
1992.13 and take t C j = 1992.65 from the extrapolation of the 
component motion. In such a scenario, A<o,i = t CJ — <fl arc , and 
we find 7j = 11.8 and 9 } = 3?2 . The corresponding Doppler 
factor is 8j = 16.4; the core-to-nozzle distance is r\ « 15 pc. 
We note that the derived values are consistent with overall kine- 
matics of 0836+710 (096). We therefore view this model as a 
likely explanation for the kinematic properties, spectral behav- 
ior and flaring activity observed in 0836+710 in 1992-1993. 



4.3. Properties of the inner jet 

The kinematic parameters determined above for the case of a 
non-accelerating (b = 0) jet with transverse (a — 1) mag- 
netic field can be further used to evaluate the physical condi- 
tions in the emitting region responsible for the observed spec- 
tral changes in 0836+710. For this purpose, we associate the 
core-to-nozzle distance, r\ = 15pc with the location of the 
emitting region at the estimated epoch of origin of B3 (< C j = 
1992.65). Then, for each spectral epoch, we find the respective 
location, 

r-n (>gf. (5) 

and size of the emitting region (using equation (3) from Marscher 
1987), 



\ (5e+a+6)/4e , Q 

■' = "'<£) t 



1/2 



(6) 
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Fig. 8. Projected location and size of the emitting feature as 
derived from the observed spectral evolution. The locations of 
the emitting feature are calculated assuming a viewing angle of 
9j = 3? 2 obtained from equating the core-to-nozzle travel time 
with the observed time-lag between the gamma-flare and the 
ejection of B3. The respective spectral epochs 1-7 are given 
next to the data points. Open square denotes the location of 
the VLBI-core of 0836+710. The increase of the size of the 
emitting region can be modeled by a conical jet with an opening 
angle 0j =2?1±0?1. 



S\ and^i are measured at t C j, which gives v\ = 17.9±1.8GHz 
and Si = 0.98 ± 0.11 Jy. For di, we use the size of the VLBI 
core, O corc w 0.1 mas, typically measured in the 22 GHz im- 
ages of 0836+710. The corresponding linear size is d\ = 0.4pc. 

The derived projected location and size of the emitting re- 
gion are shown in Figure | for all spectral epochs. The open 
square refers, in Figure |, to the back-extrapolated ejection 
epoch of B3 (associated with the passage of B3 through the 
VLBI core). At earlier epochs, the emitting region was opti- 
cally thick at frequencies v < 22 GHz, which is reflected in 
its negative location with respect to the VLBI core (implying 
that the emitting region was traveling between the jet base and 
the VLBI core). The region was expanding almost linearly; the 
resulting opening angle = 2?1 ± 0?1 is in a good agree- 
ment with the observed opening angle measured in the imme- 
diate vicinity of the VLBI core (096). During the expansion, 
the magnetic field decreases from w 0.2 G to w 0.05 G, with 
the core magnetic field, B t= \wim ~ 0-1 G. 

With the jet parameters determined, we can now discuss 
the energy balance in the compact jet of 0836+710. We adopt 
the formulation of Blandford and Konigl (1979) to estimate the 
total (kinetic +magne tic field) power of the jet, L to t = 2.0 • 
10 48 erg s _1 . The corresponding synchrotron power of the jet is 



L 



syn 



2.5-10 47 ergs" 1 . 



The data from von Linde et al. (1993) 



provide an estimate of the energy deposited in the jet by the op- 
tical flare, E ovt w 2 • 10 47 erg. The maximum observed gamma 
luminosity of 0836+710 is « 1.8 • 10 48 erg s _1 (Mukherjee et 
al. 1997), but it should be made smaller by 10 2 -10 3 to account 



for Doppler boosting of the gamma-radiation. Comparing the 
optical energy output and gamma luminosity to the jet radio 
power, we can conclude that flaring activity is not likely to 
play a significant role in forming and maintaining the jet. On 
the other hand, we find that the average isotropic luminosity 
of the optical flare is similar to the isotropic radio luminosity 
observed at the epoch of the flare (ss 4 • 10 45 erg s _1 ), suggest- 
ing the relation between the emission in the two bands. In our 
view, this supports the schemes in which relativistic electrons 
(responsible for synchrotron emission in the radio and optical 
bands) are the primary radiating particles in the jet, and external 
photons (producing the observed gamma-emission) are Comp- 
ton scattered by the jet electrons (e.g. Blandford & Levinson 
1995). Then the energy densities of the synchrotron and Comp- 
ton emission can be related through the Doppler factor, so that 
uc/u syn ~ k e 5l~ a °, (Dermer, Sturner & Schlickeiser 1997). 
We assume equipartition (fc e = 1) between the external radia- 
tion field and the magnetic field of the emitting region. Then, 
for a typical spectral index of Compton emission, ac = —1.0, 



we obtain 5 



16.4 derived in section X.2 



15.3±H, 



which is consistent with the 8\ 



5. Summary 

High frequency VLBI-observations performed shortly after a 
strong gamma/optical/radio outburst in the quasar 0836+710 
revealed the ejection of a new jet component, which is moving 
at a constant apparent speed of (3 ~ 10. Back-extrapolation of 
the motion yields a time of ejection of the component shortly 
(< 0.7 yrs) after the time of the gamma-/optical flaring. The 
monitoring of the variations in total radio spectrum showed 
a synchrotron self-absorbed spectral component propagating 
from high to low frequencies. The path of the spectral turnover 
of this component, which we identify with the ejected jet com- 
ponent, is consistent with adiabatic expansion at a slope of 
S m oc u^ 6 . The time lag between the outburst and the ob- 
served component ejection can be quantitatively explained by 
strong synchrotron self-absorption at the jet base, which pre- 
vents the detection of the radio emission from the component 
at early epochs. In this model, the jet radio-core is separated by 
« 15 pc from the jet base. The jet plasma has a bulk Lorentz 
factor 7j w 12, and moves at an angle of about 3° to the line of 
sight. The corresponding jet opening angle derived from the 
expansion of the spectral component is </>j = 2?1. The de- 
rived values are consistent with the jet kinematic parameters 
obtained from VLBI monitoring of the source (096). 

The high gamma-ray luminosity (Figure [F]) of 0836+710, 
which surpasses the IR/optical luminosity by at least one order 
of magnitude, the rapid and apparently correlated variability 
seen in the gamma-ray and the optical regimes (Figure ||), and 
the ejection of a highly superluminal (f3 app — 10) jet com- 
ponent after the quasi-simultaneous optical/gamma-ray flaring, 
strongly supports the idea that the gamma-emission originates 
in the inner parts of a highly relativistic jet and that a broad- 
band correlation - perhaps extending over the full electromag- 
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netic spectrum - between high energy outbursts and ejection of 
new jet components exists. 

The observed multi band energy output during and after the 
flare in 1992.1 allows us to conclude that the gamma and op- 
tical flaring activity are not sufficient to maintain the jet total 
power, but are connected with the observed jet radio emission. 
In this scheme, relativistic electrons are responsible for syn- 
chrotron emission in the optical and radio bands, and external 
Compton scattered photons produce the observed gamma radi- 
ation. 
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